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(jktoal Positicfpipg Svstemfi 

The present iDvention relates to developments intended , to iocrease the accuracy 
obtainabid ftom gjbbal navigation satellite sys&ms (GNSS)> 

At present tfa^ are tv^o publidy avdlable GPS systems, known as NAVSTAB, 
5 owned by the USAi and GLONASS owned by the Russian Federation, tliese have 
been in existence for aronnd tnro decades, but tn ^ near future it is hoped that the 
European regional angmentation of GPS wHI start to provide its services^ followed 
within a few years by a Buropean system under the name of OALDuBO. 

The existing systems have be«i progressively refined so that ucdng a differential phase 
10 in^lementaddn a iocaiional accuracy of less than 2 cm can potentially be achieved 
over a baseline of 1000 km, but widi a cost ia computation and in the time taken to 
determine die location. Real time or near real time measurements have a 
cone^ondin^y lower resolutioxL, and at present the requirements for hi^ precision 
mean thai additional augm^utations ore necessarily employed to supplement the 
15 GNSS informatiotL Furthermore, these could include a receiver takmg measuremeiLts 
fbom many satellites, to all diose visible to it whereby to calculate an over- 
determined position solution and rgectmg inconsistent data to improve the accuracy 
of the position solutioiL Such a system may use data fiom more than one constelMiOn 
of GNSS satellites, GPS and GLONASS. HoweveTji a significant em^r arises fium the . 
2 0 inability to model accurately the delay to the GNSS signals caused by the troposphere. 

One of the largest identified sourcra of error lies in tiie effect ttat the troposphere has 
on signals pitipagating therethrough- The troposphere introduces ray bending and a 
signal delay which is influenced by a number of meteorolo^cal fectors, but 
particularly water content, 

25 Ttaditionally, ^ delay has been handled by the use of global tropo^heric delay 
models. This weath^ data is too Isxge to be sent ov^ communicadan systeo^ that 
ara available to mobile users. 



1 



. .sur. riHNCHESTER ■ jq pptcNT OFFICE 



However, tacetmy an attempt has been made to mtroduce nmnerical weather 
prediction (^fWP) into the GPS process, whereby local estiroates of Iropospheric 
delay can be provided. The intraition was to use high resolwiion KWP fbr 
determining the liopospheric zenith delays for the estimated position of a GPS 
5 receiver, instead of udng a global delay mbdd, which ddays can subsequentiy be 
. used is the positioiiing process. 

However, even if the zeaiith delay is known, there will be a need for a user to be ^le 
to coEcect fax other directions if increased accnracy of measurement is deaied. 

The present unrentian provides a method of obtaining data Jfor nse m a satellite 
10 positioning system or GNSS comprising the steps of genstating a first set of less 
accurate troposphBric delay values firon a first model which h known per se, the first 
model bdng based on a first set of parameters, generating a second set of more 
accurate tropospheric delay values firom a second model based on meteorological 
infcnnation, and developmg a modified said set of parameters for nse with said first 
IS model so that it then provides a set of delay values substantially ^oxnnatii^g the 
second set. and expressmg at least one of said modified set of parameters as a 
coirection paramet^ for commimication to a mec. 

M an embodimejit of the mvention, fiiese parameters are difiersmces with those 
obtained firom tfw fifst inodel and are transmitted to the receiver. 

20 In one enibodiment the finrt model is based on one or more non-meteorological 
paramater(B). These pazametei<s) comprise at least one of (A) time of year. (B> 
latitude and (Q attitude to generate an estimate of the zenifti tropospheric delay. 

As particularly described, the non-meteorological parameters a» employed for 
detennming a mapping function relating tropospheric delay at a given elevation angle 
2S ID the zenith delay ganecated this first model. A m^ing function using a 
canthined fiaclion function hased on ft.e three nonrmeteorolo^cal pazameters such as 
used hy Niell may bo used fbrtUs process of geosiatbg &e mapping fimction. 
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The second model may be based on NWP or real time meteorological data, wiuch is 
used to calcnlate the zemth fiqpospheric delay. Currently fiie cell size of NWP data is 
of the order of 12 km, but this is expected to xedizce in the foture allowmg finer 
detennination of the model and mapping fimction. Using the NWP or other 
5 meteorological data (collectively NWPO) a ray tracing techni<iue is used to determine 
the path of the satelh'te signal through the troposphere to the user and hence estimate 
the delay from a direct path. Knowing the delay, a second set of mapping function 
parameters may be derived by fitting the results to the format of the mapping Amotion 
used in the firisit model. 

10 Rjecogoising ttiat the two models have generated two values for the s^mie set of 
parameters, it is only necessary to communicate to the user the differences in the form 
of corrections (ag, as actual numerical difierences, or as a multiplier of the value in 
the first model) to the four parameters estimated (zenith delay and the three 
parameteis used to generate the mapping fimction) in the first models although it 

15 ^ould be possible to cammmdcate the parameters generated by use of die seccmd 
modeL 

Depending on the nonrhomogeneily of the troposphere it may onfy be necessary to 
generate a subset of She parameters althou^ the most important is the zenith delay 
followed in order by A, B and C. 

20 To be of USB the datamui^ be commrniic^ted to the user within the time far which the 
NWPO is valid. As noted above the NWP data is defined ovct a geographic grid. If 
tiie standard deviation fox titie change m z^ith delay for a cell in the grid in a period t 
equals the inaccuracy in the second set, then for a workable systom (where the 
me$^surement period is less than changes in troposphedc conditions leading to 

25 misleading r^ults) the data reduction n^ds to be effective to reduce the data to an 
amount which can be transmitted over the capacity of the data link in a time 
subsy^tially less than t, and preferably no more than t/4. 

The derived correction paramt^ers are distributed according to location over at least 
part of the earth's surface and may be communicated to at least one geostationary or 
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non-^geostationaxy satellite or tectestdally based systeiiL The GNSS satellites may 
also be used for this coxnmimleatiotis link. To liable the use of low capadiy data 
links it is piefored to compress &e distiibuted data prior to transmisdoiu Such 
reductLOBS may conveniently be perfonned using lossy knage canQmessioxL (such as 
JPEG) or sub-sampIiDg tedvucLues. 

Furthex details a^*^ advantages of the invention will be evident upon a reading of the 
appended claims^ and upon a leading of (he following article. 
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ABSTRACT 

Witb the nsKiennsetiaa tsf (^S ^od the nstednc&in of 
Galileo, icnpo^li^iic enors are ' to become 
ioeieasiiig^ s^aficanl piopoxtiQii of &e GNSS ttscr^s 
ciror bndgci. Thb p^r invwdgaies chc spadal a&d 
tenyxual propenies of zDcteorolosi^l Teiitocs oa ttis 
tropospberic tma cxpcnaxcd hy GN5$ weis. A 
KgiodBl ra]-D*m cropo^hed^ conmioa service coold 
be provided fiom rc&siQid^ profiles, deiived fiom a 
nunusiJal wcaOer ptediedoii aiodel. Aa aaaljpsls is 
described of Ibo bandwidth reqavcxncnts for Tcaktiioo 
troposphnic concction dbsamsadaa and of the icadc- 
off between bandwadlh and accumcy. 

Whilst the dry atzoo^heit vteijib£Iit>\ botb spatially and 
tenp>r3l!y is stualL water vapour variabib'tjr can bo 
Isse. High-icsoiutioA mcso^Cdlc Tinmi-^fT yi wathcr 

fields cittblo aiBl^^ of the abfupt spaiial Atid tfiopord 
cbaages ofa rmoitphrri g zsoisture associated with weaihor 

ray-tfS6ut^ techznqiiip has 

been ^vtHopcd to coK4>i»ft nopo^twtnc delays fiom 
mesosealc pumncal wither pppdictmi &lds for 
defined ismmsb flievatloo and oziaiatb aisles. 
Quazitified cxan^Ies of ritfOpl ^tJal, ioeludio^ 
azfimiTfiBl, «d]rnimetxy pud (tsnipwd t3COpoq?hecIc delay 
vaiiabflity ate given and sutisttaa aadysla of tfaeir 
likelihood presented 

Once the tufpoaphsrio delay can be established, die &sue 
f jbtz caaoprtSsSm of regipnal fioppspbcric 
in tewMfa B to torn a prtetieaJ Mxvfee needs to be 
adcbessed Diffmsct stuat^es for disscmlaatioii of 
rcgionad troposphcric oonvcdons are discussed aad fhp 
tadv^oSf betmea scciffasy And coocctibo service 
bandwidth described. 



1. INTRODUCTION 

Tnsposphcric delay CdJKCtfams catt be derived from a 
jHion atmospheric xnodcl, an atnio^pbaric model ufilkiiiB 
xcd-hmc xneteorolosical data Of* {ivea dufiScient txin^ 
estiTOied from GNSS code and eaxxMf observ^tioas. 



lias paper investigates (he pofemlal of Nuffleifoal 
Weatisr JRredlecioa NWp. xnodels 19 pnmde 
fropo^hcdcconectioiis. / 

The 3«D atmospheric pxofile mfonDaCio!! resultii^ fiom 
the pn?ce5s of numerical weather piedic?tiOa oSbts (be 
pomtial of pnnodiqg a les^pqal or ^obal Iraposidieric 
wixecdoa service to GNSS usees, la eAbct aa accanite 
atmospheric refiactznty field, and dicrelbre tropospheric 
delay. Is sbwply a by^rodort of (be owacripal weather 
prediction process^ AD accurate estiiaaTe of oopo^dieric 
dehiy can be obtained by Toy-trasiftg thn^ti^ fite 
rcfiactiYity field, the accuracy of a NWP-derive4 
tropospheric correction is almost entirely d^eodeot iq»Dn 
dw absUQr CO deiemrino BUnos^faeric tvater yvpmr 
content 



2. SDRAYTEUCING 
2A Pr^UieatloD Theory 

The speed of pcopagatioa of aa clcctroinagxietac wyver 
throujE^ 9 mcdiim C9D be expressed in terms of cbe 
re&aedve index, n, defined to be die ratio of the speed of 
light thiooghfiee-space la tiMt speed through the medium 



c 

V 



CM) 



where: 



n Mtibc g e&aetSve index 

cIs tfke qpeed of light jn flee spaec 

V is thepropag^tite vdociiy 

Tbs GNSS tropospbcria time delay, ignOxiiig lelativistic 
a'fiscts^ kf dc&Eed to be die p r opsgatioa im^ of the 
C^^SSsagnal fiym die sateSite to the user mintts die fiee 
i^iBcOpEOfij^gatiOo &i3es 



Vmr Iter 
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X is tlie dlctanee bIoo^ the impa^tLDn path. 
The fixst integrol is aloais curved 
propa^tton patb; tb9 $ecQo4 infcegral is along s 



The dif&rential equation dwaibing the cttfved ray path 
can be ejqimsed, in cwesian coonfioaiest as PJ- 



Vn 



r /&) 15 r&8 vector describing the ray p^itli. 
jr to the Icxt^ of iBe evivcd rsy pBtii up tff 
A b the lefiactiviiy sCftter iitld, 
7tt« s vector ficH is Ite grad&eiir of 



dif&csirtial oqu^tlon era be esspaadbd as 



IF 



A first order csdlnuy difiRxctttiftl equ&iioa wiiiti 
kno^ initial values ean bo solved usuj£ imrksedca] 
methods: for example RuDgc-Kutts or Adaij^-MouliOD 
methods. Higher oider difiEbreBtxal equatiDSis ew be 
solved Qurocrieiilly by renvziliiKg thtox as ao equivalent 
system of first order equatfoias. Usizig the mb^ttttion/i 
= r and » r* (ihc {icst derivative), the ray path 



diffexential equation (2-4) cm 
cquivzdcnt sysl^ of two ftst 

. dt 



be esQiijre&sed a$ aa 
diffitf ential ^■i*^f * nt 



The dengnfninbg of the my path tbezefoxe ax&oufits to 
the sohitlon of a system of two OD£s with initzal values, 
Matlab'sfuneticaede45 |]$«d ta $oive ihe piDblem; 
A Rung^-kotfai n»elhod wid ad^^ve stq;i eooatzbl 
coisdsteitfwifiiUfierdeflned |oIeiaooe& 

The mimical weather psedtedoQ liclds to he used are 
c^qxressed in a ^pbeticitl coordinate ihUD9» it is 
cM|7at»ieinA% cozxVtAieBK. dmefm, to generate the 
lofivctivi^ gtedteat in sphencol epotdraates ^,£lo)« 
which can be cuuvcjitgd . fiato local etnveilsDBar 
coordinates (ii»v,w} using the followiitg iiaasft watioto* 



A fiiither £0i4tianat fnmsfornttdoii is then applied to die 
local cniv«IBiker coofdbiate Ihune to jgive die gradient in 
1 fpc*d cartwiro frairw (a.y^ ia Fig. «i£table for 
nontericalty solving equs&tion^ 2-5 eMl 2*fi- 




H^Z^l RcfiaciivixyGeoxiaetty 



The tropospbeiic delay can be compoted as: 



(2-8) 



where a, b and c ace M shown in Tig. 2^2^ Point b 
CQizespoods to l^e pouvt at wfaicb xay ctnvature 
. sad ic&actmty is negSgible, hi lliis work 
aSStttoed to be ^ve an altitude of TOlon. 




Eg. 2-2 Ray?3itbCeotaneay 



EmuplePMpaeatloDF^Ilis 

An example prppo^ticnL pafih f<ir a saaeditB at an 
etevaiioa of 5^ is sbown ia Fig^ 2*!^ The co*ordiafites 
used arc and >^ as shown la Flg^ 2-2. It is soted that 
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•pjodty of tbe cmvstun; of the ray^pfttli occurs lOCK- 
InftonifhEUscr. Itiftlheroutionoflhereihctivit/ 
rnt vec cor Cdd, caused by eutb CQxvature« «4»efa is 
the ddminant «»u$e of my curvafucCi tfttbtt ibaii tb& zate 
of daqgc of rc^ijncdvicy scalar field. Jhc sppxcent 
elcvadoa angl^ the propftjS^on p^dli wgle at ihe user 
locaiiofl. is &ppraXiil»Cdy 5.19^ ftr » true elevatkoa 
anglw fiom SPteUite to iser of 5^ 

Uis shaw^» in FIe^ 2-3, fbat an eadicnie sonosplicre CTBC 
slOO) doc$ noi xesuU fa 9 ^^Ficaai claii($c in s»y ftngle 
entry into iiie trpposjpbm ^boerore^ tba assus^scioiiof 
fbcaddStifia of £rO]MS;phmxo ionospli^ delaiTS is valid 
AS i3ie couplisag dfcxt b sniiAml 



o 




o 



Fig. 2-3 EKffi9)IerSAyI7«v|ati^oaal3^ 



Example dddys and ^i^idanal p&Ui length due to niy 
cnrraturc are given in f2£. 2-i4. Tbe incxease in paiL 
Icogili dnc t9 i?iizvatctt is apptoxinaicly I9cw Ibf a 
satellite at 5* elevatim. 




2.3^ Equivalent Angle Ajiprou'mailoa 

The roetbod of deternuou^ the. tzoposphex^ delay» by 
numcifcaliy solviiiS a syvlem of ordinal^ dSflbteoiial 
(quatiacs^ is coJiP^uisdaBslIy t3q>ensiw. A bstcr 
method is lo inteigiale the refiactivc mdex ftloi^ a 
$trai£bt geometric path; however this iiite^atlou wHl 
result hi as oveire^tlDSation of tiopo^hcric delay; la 
reslity tbe ray betuSs so as to ivtlrtifahc tiie propwitvon 
tnrv (Fetxna's piiaciple)- 

There fcxisis an eqinvalem elovatioai fti^e, tying berwcen 
the true (geometric) and sppatt:^! elevation angle, tbat 
results in tibo straight padt int^gxa) being equal to &he 
!nte.gisal along the actial corvcd path fiom SAieOlxe to 
user. Hit elevuian eorreetion to bo added to the uue 
elevazioR to obtain ibc c^aiv^dcm vasies with user 
sft^tudc; but IS largely unchanged '^fh aimito^lbtn. (Ftg. 
2*5)«. The standard devution beiwcea the ct^azvaleol 
sualght line integral and the lay-lxaced troppspbcric 
for elevations between 3 and 10 was fimvJ te'lir 
appioximainiy t).049&» 



0.12 
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9 m * ft 10 12 M 



Rg.2-4 ]$KM»ple D^:^ and Excess Q^TvedPaA 



Fig. 2^5 Cdrreetion CO ItaeElevadon Angle. 

a. rnJME!aCM#WEATH£RnEU>F^ 
3.1 Istradnetioii 

NttoaetioBl Weather PrcdictiQa xnodels ftxrecast Ihc 
cvoimion ofaui^^kedc physical processes by applymg 
^ovemii^ eqpatims^ including the conservation of naass^ 
mMMBAiDl and eneisy* three-dimensioual SUlds tit 
coninuieus variables fnflnrWng huzzudityt joesstir^ 
len^petamre and velbciry are ammeri^lly pfoc^sed and 
ae^coiologieal fiatnres, reh«Jing wcaihCT fioots, arc 
secondary deriVod pwperties, A iraiiety of 
iY >f>^^FwsggM?Ti>g i-»rw ^ ynpnt Inio fbs smpcticdl oiodel 
laclodlng snr&cc^ MJiosonde and 99lenite ofasetvatiDas. 
The water cydc is fxiedeHed Sncfaiding the eRbcts of 
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xnoi^nifep sea sisfiice tenapBiatixrc. elaud 
And pxccipilAliotL 



single Cu^tenn eTcpresmo, s^*^^*^ ^ C^l* ^ * 
refiaetivily accural of (X5%: 



o 



Tl^fft m & variety of co-osdisuie sysEcois in use for 
Dnmericftl wearbcf |>rediccioQ; they can be either based 
upon 1 grid of fixed poifiis or spectral w«»ffid«at$. 
PbysicftI hcigbt is mcty used as tbe venic»l co-ordml? 
ui NWF: picfisure, sigBfti-picssurt and hybnd co- 
ordizL&tes arc tooxt codvcoxcol SigouL-pjrcsscire vcrticdl 
co-ordinates define a to at surface of csoni^l 
pressiuK divided by sutfxce pressure find tims haive a 
lenain following property. Hybrid co^oxdicsites artr 
tenain fbUowiag at Jow tevcISp but teod cofltioqawsly to 
pressiBo $url^c^ 3-t upper Isveb. 

Nauerical m£>dds ca^ be globsil or of Limited area. 
Ximised high^-resolutioo psodel^ are ofieo termed 
mesoiscale models as tbcry reflect mesoscale 
i i i e tcQru i<>gjeal features: weedier pactexns of less tban 
lOOkoiinsiSC^ 

Tlie nsesoMale data used i± ibSs study was abtftiaed firoxn 
tbe British Atnuspbcilc DaU Centre wrdi permission 
fiom (te trinted KiASdero Meteorological Offiea. 
UKMO, The UKMO iiasosc^ Unified Model Mse» « 
hybrid vcrii^Al ccMrdioate system with a spatial 
tesahttioa of 12x12 koi. Tbe model covecs the xtcA 
' shown hi Fig. 4-2. Specifie bnibidliy sa noodeUed ax 35 
hybrid levete; tenq^enture is laodelled To v^y- 

tzaee duougb die amocrieal westhcr Held ite hybrid 
levels ar? inverted inCo physiOftl bcighx by assuxoSng 
hydrosudc cquiUbriuml The zfCipaet of dcviatiozu firom 
iQfdrostatic equilibrium csiased by mouncalnovs terrain 
on lunar laser RUigc msasureniBAts has been estilTttted to 
be less than 1cm for etevatton ffiagles &s low as 2(r IS}» 
the inspaCT due to oonveetive Sloniis bt sot knowo. In 
assuming hydzOSUtie eq|uilibriuzn« local vahies &r 
gravity are used to reflect £;ravita.tionaI vaiiatioo widi 

Ooce ihe NW7 fields are fixed in 3-D space^ the dat^ 
W rc-sempledr USiOS a romtwnntten of Xfoear and lotg* 
linear hiteipd!adon» at va^ pl^'cal height, Ihe 
atnaospbere -was cxlrapoldttd from ibe top of the N Wj^ 
field (appropdmately SSlaco} 10 an alcinide of TOScm 
n^nnmin^ a log-1in«r lebtionshfp bctwi^ i«£^ia4viiy 
aadoltitwk* ThcatmdiQifaeitteaath&nberqpkresenK 
^relbcttvliy fi^li^ fiont i^di ihe refineiivicy gfa^tiM is 
numerically coinpntedL 

33 Atmospheric Refr«ctfvftar 

The atmospheric it&activity can bt di^ved j^m 
ataoo^p^hsric preasuzep water vapOcir partial pressure and 
tcxspaatuie fietds. Atnxxspberie xtsScktiMxy cm be 
divided into dxy (bydrosta^) u»d ^vet con^nents. A 



where: 



A^=10H«-0 



(3-1) 



/iTis the leAsCcivity. 
P is atnaospherie pleasure (millibar) 
£ is the water vapour pressure (niilllb&<} 
Tis tbaapmami^ (Kelvin) 

A r^vi^w of sophisticated thxec-tezm essp^essions 
including gas cosx^nrssabilitr ef&cis. possessing ptmt 
accuracy, is given in 

3.4 PwtureNDineriealWefttherPr^ictfon 

The rapid rniproveine&tt la cQispotaGonai ptinvtr wod 
advances in nufflerieal tecfaniqiues has seen the ECMWF 
BX}dcl*s horizontal spadal resolution sedufie fiom IVOkm 
to 2ti1an in (he period 1979-2001. This process is 
pepected to continue a»d it is predicted ihat by 2015 
Umxxed-ma model rcfolotion could be 1km and tbc 
global nsodd leas than IQkm [I I]. 



4. MSSO^DAIJSMETKOAOU^lCALEEATDaES 
4.1 Relevant Mtieorologie^l Features 

Meieomlopoal features tbac possess a large apalial 
juid/or teix]poral variation in trttpo^heric delay wiJl 
Inspzci itie accuracy of NWP-dcrived txopospherie 
CQiTcctions and fb? bandwidth required for 
dissenuiiation. The ten^roi?] and ^atial variation fn 
hydrostatic refiactivity is jgenezally stnall, wbtfrcas 
oeteorolo^al features associated with rapid rhanga in 
atEoOiq^erie nwismre sig^iificantly Inapaoc the 
Bccaraey/bandwidth relatioxisbtp. 

ITic JjrtJjal and tenipozal scaled of cenaia weather 
sysiexm £tfe described in Fig;. 4-1 (Crom [1])r 
MoteOTolosical features smaller than the resolution of the 
oumerical prediction oaodel wilt not be aceurately 
reflected in itie NWP-dfitived vcqpospbeA: correctfon. 
the 12k}2 km tesoludoa fklds nns able lo reflect 
characteristKs of we?ither Itpnts^ however smaller 
feanires such as iMal stoms and comadoeS wobld not be 
well represented spatially or tatnporally. The extent of 
Utt weadier fcsitures in altiiude is critical in assessing 
tbeir rdcvaaee to trqpospberic dday asoroalies. 
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]^4-I Spatial Sc^le ?nd Lifi> Span of 
Mcteoiologicd Systems [1 j 



4J1 W^AiherAronls 

A weaiber fioni maAs the intni%iC9 Mwe^ air masses: 
defined as A luge bod/ of sir whose pb;^lca1 properties 
m Isrgdy unifbnn hoiizontalSy for liimdrcds of 
Uoroeires. Tbt &oax can inajk tbc otrojjrco^e ofibtrupi 
changes in Atmospbenc Tnoi^tiue, ten^:catur« and 
therefore r^^ctivity. Fronts vm bv divfJed inrc? thice 
classiiicatiozis: waraii cold and mlud^dr 

Ezth frontal zozie is uroque howevisr ibcy do posses 
some geneial cfaaracteristitis. Fxoutal zozu» tend to have 
a width of a few lOkms. A -vnam from repr^s^nu the 
leading edge of An adv&aemg wann sir mass, axtd 
gcxMnaUy has a gentle s-lgpe (ix)i;1ina^oilX ^pcally Qjr - 
1^. Cold irojaxs, th& leadic^g edgea of advanduag cooler 
atr masses, are xypkaXiy indizied at 2^ Ano&ar 
cfaatactcmtic of the Boat is ix^ relative speed, cotd fiml^ 
mi^b |a$icr than waxm fioms. iip to 50kzate 
cojx^raicd to ZOkmilur respectively. Hence, when, a cold 
from catches with a warm fiont ft radeiculs the 
wanner less dease air and fba cesulfanu mene eon^lex^ 
featvra u temwd an oceluded front. Occluded fimis 
contahi elmactariacies that belong to both cold and wann 
items. 

Use roqyl r^pfd change m trppospheric delay is likdy to 
occur when, saidlite eleradon and iiom iT^itrif>f^'nn ^zo 
equal* Goiexally^ iii the ll(fC, jh)mal system 
to 50 kgometnes per hour aod cap TWrtt bi aniHh delay 
variadoDs of ScroOiDiir [ 12}. 




Fi8.4>2 Synoptfe Chart WeadwrProtit Feature 




Fi;|.4-3 Meteosatlnaso 19/01/2002 



Aa cxanaple of a weather Eccaxal aysleci is shown Ut the 
^■ttoptie chart bi fig. 4-2 (obtauicd iteoi the 
Wetteizentra]e website)* Tin: n^usn witUa die 
bighlffibted Uidlled Modd twmgi area, hndod^ 
pcchxtcd, eoid, varm and elevated eold fioals. The 
fionial sysiesi is associated wid> abn^ changes in 
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Fig. 5-1 AsinmUhal Asynuneiry}^earcd£ii^c 



atmospheric moisture comm% whiA be seen as aa 
abnipi cliange ia cloud coverage, Ffg. 4«3, a Mctcosaf 
image obiaiaed iiam Dimdee UAiversity; aad specliic 
humidhy. Fig. a NWP-mesoseate im8$e. 



5. AZTMIJTHM. ASYMMETRY 

The scc:Ur^y of tbe widespread ftSStunpUon thai ihe 
aoyM^slb»c is syu uu utiic in JtzimuKh can be tested by 
^xaimwg the vflzi&don ia troposphecic delay for b given 
elevation BBgle. A dati $et. eonsprising UXMO 
mtsoscaic NWP fields in iwo^urly isieiYAls fct ihe 1" 
and JS^ of each moatlx over s period of quc year, was 
ised t9 generate assimutlul ssymmetiy anplhude 
statistics. Fovr uset loeatiozis were used near the 
mcsoscale region ventre: each locatiofi -was tested in 
uscrcmimts of 10 degrets of aximuQi at an altirudc of 
500di above meft-sea-levcL Fig. 5^1 show» the 
percemsgc vuiaticin inlnitpwipZkcric delay fiom die izicaii 
delay ibr dovBtiDn Wgles of 2> S» tft and 20**. The 
BisBkdaid deviation of azixnufbal asymmeiry express a$ 
a pvrcenjj^^ for eacb devBtloa is aZ]« 0.07 find 
0.03% respectively. The amoimt of aslmuthsl 
asymmctzy is fiznall (q relation to the current accuracy of 
die NWP-derived trOpf>5p3ieiic coiri^oiis: a standard 
dcviadon of 0.4%. Por usesfi requoiag d)e hi^liest degree 
of aceuracy, BzixnUOial asynunecry infbzxnation -would be 
of value^ however lor imst vsbis, the cost of 
disseminating T^rinnifhgt asynunctry Xutfornttuoii 
Quovcighs the £znaU gain ioaecu^acy. A« the accuracy of 
NWP-ilenved ccnnctions improve^, maimly due to 



improved model resolutiaa, azionitbal aayxmsetry is 
Ukdy to becGim a more sxgnlGcaot Unu'ting ihctor. 




2000 



-I .^M M JOA ^3 n 0^ O.C> €g 




10" 



■1 



-Oa 0 OJf i 

20^ 



-ej -9,6 'Qii -9.3 o o»a 0^ 0^ as i 
Aslmtilhai Afiymrnat^ Parcanfeaga 

Fig. 5-2 Azisftuth&l Asyxnmctry Hiatogram 



The distrftution of the measure of asymmctzy has long 
tails. Exiretae leveh of a^ymnaeiry occur when a 
wither fioni is oveirbead Fi$. 5-1 consists of a specific 
hmnidSCy honzontsl slice at an altitude of 2kin, & tyrsdRc 
humidity vcztical slice and polar plot indle&liog the 
percentage of azlmnthal asymrneiry. The polar pic* 
displays the percentage change in dday IStom the meaA 
delay for i(bat elcvsiion. The conccntzic zuigs are for 
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•ioo iwigto of 2, 2, 3, 4, S» 7, 10 aod Tlie 
of suyrninetiy causal by simple weatlier fronts 
0 accuxaidy modeUed » a simuoid. Sinusoid 
paranu»ci5 could ihan used by usm requiring the highest 
accuncy gf troposphcdc COitectians. as suggest m (7J. 

A ELEVATION MAPPING 

in jnoddling trojpc^ph^c delay it i$ cpuvcrfcnt la wJaic 
nopospbcrb d«ldy a given elevaiioB angje lo ihe 
ztmUx dday (c^} °i^>>n^ unppm; funcdao 



(J 



Tte hydtostatic seniUi dday c&n be aecurastly iradc^lcd 
givea a surface preswe ZDcasiuwQait bowcver ibc 
WBC aea!ib delay cftn dOC be Mcsurately detefmined from 
sor&ce humidity n^surezneuts, 3S iboy are gooi 
rcpfcscniadv^df the Above atrooapihctB. Tb&expres»m 
for the tropospfaerie deby at ft given dcvadDn asglt eaa 
be defined as; 



It is fioted that eaxe mwt Uken when applying (he 
fiimp]ifictitiDxi of rb? supecpasiiioiiof bydzoatatic and wet. 
atioosj^iecie delays: As propagaliOn path b deprndeat 
iqsonbodi hydrostatic aikd wet compodftA(S. 

The variation of tropospheric delay with clcvaiioo angle 
(Pig. 2-4) can be efficiently modelled by a contimied 
fiaction eaqpannon. NiaU l^drostalic and ivet icapping 
functiotts [10] p0StaS« a bi^ decree of a^cwacy wilbo\it 
ifae need Ibr priot ineteocologleal {nfaiznation. Nidi 
tuning iuBcdoBs use a cDiuintsed fzaciioa Itiscdoa 
based tn three paiaxneiers a, b and ^ which are 
fiincrion of tttw of year» latitude and altitudeL Niell 
tDftpping fitrrnKnw^ are lOcely to result in sigfiiflCMK 
addBtional eirofS coirpaied to the aeeinaey of NWP- 
derivtd zenith coneclioini* A mckie accurate di^eci 
mapping faas^xim can be used baaed upm NWP my- 
tiaced resales [4]. 

For die reasons of disacxmnatios) e9ic{eDcy» it was 
decided to m a duee-ienn condnucd fraction 
approasxoation to tba totaS (wet sod dry) mapping 
function deiemiified fiom xay tracing NWP fields 
(equation 6*3}, P3nxoetenao.'ABaxidco9untanCfoaiia 
pritm HkQ mapping foactioA for die total {hydxosiaiSc 
and wet) troposphede delayi uid tie a fimctionof time of 
ycar» ladmde and aldtude. The p&tamectfs ^« 6&, end 
Ac are a totrtetion toim ta be broadcast. The correction 
pazametcrs are detemEned fiom -a I&fing prioces& audi 
tbai die- mm of squares of the restduals between 
cquatioa (6-3) and the ray-Mi^ ttmdO is rB b dnm a ed . 



&*ch an tipprcffdnndon was tcrond to be in agreement 
widi die ray-iiaced results (typically «0.1% rwidual 
euor). 



7. S^ATlAI, ANBTEMPOKAI- VARIABrLlTy 

The tsiie of information required to ftCCUnRrly 
disseminate faroposphecic correction infonnatK»El« OA ^ 

ceginnal basis* is dependent Q|>on the spatiat and 

i^n^oral variation of tropospheric delays. 

7.1 Temporal Resolution 

Th« Sn^pact of temporal resolution on the accuracy of 
KWP«derived tropospheric corrcetions can be analysed 
by cOd^aring the change in zenith dday at time to and 
tQ^^6L Vsing the data act as vsed fat the asinmBial 
asymmeny analysiSi the changes in zcnidi deliy wete 
conaputed for At vohies of 1, 2i 4 and 8 how tor ^ 
points in the NWP coverage ajrea- The standard 
deviation of tempoiaJ vari;Ktiont expressed as a 
pereeniagc for each &I value, is 0.47$, 0,594. 0.79S and 
1.065% resprctivcly. TI>e cOTTWOwSng Wstogrnms ato 
given in Fig* 7-1. 



e 




Ibour 




2efM Delay CtmnQC Percentage 
Fig. 7«l Ternpoini ftesolution HistOgiram 



The in^cc of teoiporfli! lesohnlon Js sigtd&cant in 
EClatien to the cmTem accuracy of the NWP-dcnved 
txopospberie cacrectionss a standard deviation of 0.4%. 
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Flfr 7*2 Front Progrcisioa 



o 



A significant error additiosal error is tbm&i^ pxoA^ced 
if teii9>oi3l resolution exceeds 1 faaur, 

TLe duwlnuion of the ten^cftl Vftdacion I»s loi^ tail£« 
Fig. 7-^2 displays fhc pticent9i£9 change In zemth delay, 
Sbr vftltids of 2, 2, and 4 Iiociih» bs a weaUifir from 
traverses the UK- The weather front has nsullcd in a 
maximum me of cfancge of sexiiih delay of 4cttVhottr 
rr^casured betv^een to and to -M hour. 



7^ ^jtthd Redolvttoa 

The impact of spatial rosolufign on the accuracy of 
NWP^derived tropo^heric corrections can be analysed 
by comparing the c}»nBe in zenith delay for locadans 
separated by fixed dislnnC6A. Using iht: data set as used 
for ifae aznnuthal ^^ymmetiy anajysiv the changes in 
zemxh delay wero coinpVTed for spmA separations of L 
2, 4 and 8 mesoscalff grid points in Ncnft and East 
dkecixons (equivalent to 12, 24^ 48 and 96km 
£6paxaticn5). The sidiKLstid delations of Spaixd 
vatiatioa^ e^qxressed as a pcrccnla^e Ibr cash 2^cial 
sepamtSon vulttc* are 0.138p 0^13, 0352 and OJSl&i 
respectively. Tl?? cmrespondspe histograms m given in 
Fig. 7-2. Sp&tialttsoIuiiOBSbJflcser than the resolution of 
ihe mesoscale znodrl. 12 Ian., are Iherefore found to result 
in sigmficant additioii trppo^beiic errors. The spatial 
voriaiion m TTopo^eric zenith deia^ has a close 
reUilonshxp to ten^r^i varintioiL As a ^veachef fi'oxii 
passes QveTv for example at '40 km/hOuT in the east 
direction, the difierence is zenith delay obsexvod at I9 
and tori-l hour is siDal» lo observed between two 
pQittfassparated by 4Qlanin the east^vest direction. 




T >0.B .ftfl ^Cij -QZ O P3 fl.4 0,S 0.8 T 




.1 ^» M 4»4 -ftZ U M as OLB 1 

Zenim Dalay Cnanga Paraentage 
Fig* 7-3 Spatial Resolution Hismgran 



8. BANDWDDTH REQUIREMENTS 
8JL Quantlmfan 

To achieve cfScicnc tr6p0SpheiiC correcliott 
dissemiM^ton* it was decided not to model the 
dimDspheire, for exAn3)>le with presstve, fiunidity and 
i^ae nte rofocaialionf but tP supply ths minimnm 
hdbnnaliDn needed for secwnlc mSlif deiemihisioit. 
Simiterty. for reasons of bandwidth cflicieney, zenith 
dday and rmppuijg; fimction psrametas can be expressed 
as coneCtlOQS to a priori models. The zeoilh dday, 
dtssermnated as a perce&ta^ coirection 10 the RTCA 
model [9]. can be coded using S hits srpanaiag the range 
.10% to ^lOA (consbienr with (13]) ivzth a tesolutipn of 
0.0781% and tbanelhre mxnimat qyaoeistion noise. 
Simlailyr mappsiag flmction eonumied fiactiAB 
inibimatiu i i can be dissexnhiaied as a conecdtw to a 
Kfell-like total delay mapping fwctipn. It is jnpposed 
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Fig. S-l Zenith Delay Izoaae Couppmsim 



ibat mapping fioiction convctilon panraTneters can be 
adequately desciibcd by 16 Bits (7, 5 sad 4 for 
puataoaezs ito, Afr aod Ac; Equatloa 6-3). 



I Sub-wnpKng I I 2S I 0.<ai% I 0.5% | 

Table 8-1 ComppeTncn Properties McsoscbIc , 



8.2 Jmfi^ CMipr<MlDn 

inv««t?£»te ihe impact of quantisation and 
coii^resBiOD oa the accuracy of dilemma ted 
iroposphehc delay correctioasr an example zeoiib dcJay 
ruip derived from UKMO mesoscale NWP fields, with 
an acrive front, was qnaatiSMl azkd CQznpressed. Ibe 
accuracy ihs recovered (disecnhsattd} GOtrection was 
Aen compared to the original values; dhr^ved from ray« 
traced jBiunericai weaifao' prediction fields. Tb£ xesuliam 
gcey^Snile bitindp Im^^e^ 24^x19^$ Bits, and three 
JPEG images, -with succcsslvefy incxcaSXD^ Cooprcssfon 
rarioA. are shown in the frrst four images of Fig, 8-1. The 
sc£e of the 34>ii bitimp ijosags; a z^ped (loss-less 
con^ressloo) bifoup, asd three JPEG iniara are dvwa 
inTaUeg.!. 



The hizmap (BluIP) aod appcd bitmap resides id ft 
quantnatiori noi^c of 0.23%. The JPEG50 image has 
con^rcssed the bitmap image fay i fkttot oi24,6, hut die 
idditioaal quantisation and compression noise is smalU 
l^s than Oj%, conptrcd to the corrent accuxacy of 
NW?-denved corrcctioa (0-4%>. It is noted hcwevcr that 
sioipte under-san^Iws every iifUl pixel, tbm using 
biliQciir iatezpolation between sxissing pixels, rtsulis jh a 
coAQiressicfl factor ot2$ and less nols^ thaa J?£O50: ibr 
this image JPEG inv&se coinpression cTTcn no advantage 
tivtt a reduced aamplixig rate. Funhe^r oompncs^on 
resulted in substantial reduction in perfbrnianco as can 
been MM 2n Pig. S-l. 



Type 


Size 
Byte* 


Co nip.. 
Ratio 




Worsi 
Value 


BMP 


141478 


1 


0.025% 


a07% 


.Zipped BMP 








0.07% 


IPEG90 


20195 


7 




1.2% 


/PEO50 


8437 


16.7 


0M% 


1.8% 


5uhnsanH>iiD8 




4 


0J7% 


3.3% 



'NotA: AoiM stuidani deviau'qa is B».*tis2iieT 



Table B-2 Coipsnssfon Prqpeities Global 



Typ« 


SiTc 
Bytes 


Compu 
Ratio 


Noise 
(std) 


Worst 
Value 


._ BMP 


2S014 


I 


0.023% 


0.05% 


Zipped BMP 


11212 


2.5 


0.023% 


0.05% 


;PBG90 


231S 


12.1 


0.042% 


0.3% 


JPEG50 


1139 


24.6 


D.093% 


0.5% 


JPEGID 


734 


3«-2 


0.28% 


2.5% 






4 


0^035% 




l-i^i-S-f-^"''jMl 




16 


0.063% 


0.5% 
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JPig. Exas^lv ZenitH Detey Illative to RTCA Modfl (%} using UKMQ GSobsl Model 



An eTwnplc VKMO glgW NWP grayscale blcmap 
image and two JPEG isQ&geSt wjtb successively 
tDcmsiKig compicssion ratios^ are »)»wn in the last tlnAC 
imsg^ of F(g. 6*1. Tlu: origtoal s^enith delay, expmsed 
8S a pexcentasp cbwgc from fbc RTCA model, is sbown 
ia Fig. S-Z. 7l)e global niedd has a resolutfon <rr5/6 and 
5/? degrees in longitude And laliludc fcspccUVcly: 
pioducii^ & global image of 432x325^S Bite. The 
JPEG90 image coxjopjcsscs the bitmap image 4 ^ctor 
of 7 (Table 8*2) with coznpresElon noise of k»« ibaik 
0.1%; tbe JPECSO ima^B results in signilibaal sddiHoaai 
zcmih dcli^y errors as can be seen In Fig, 8-1 . For the 
global image cfae sub-sainping technique, with a 
fci»ffl pf»»sSon factor of 4, zcsiills m fiubstaatial addittbaal 
compmslm ^rxofs- 

AlsuxtlxBg that mappixij^ flmetioA pamxieter coirecnotis 
can be siznflarly compress 

0«-* Bandwidth tf.g. 200 J9il/S global dAta dtssezEunated 
in lOmnute^ 

9. SUiVlAlARY CONCLUSIONS 

An accumte atmospheric relractivity profile, and 
tbattSm ttctptrsphaic delay, is singly a by-product of 
die numerical weaifaer piediclion process. A vector 
C^cuhn? bwd r^pproach reduces Ibe problcca of 
propsfiatioa pa^ dermnlpadan, through the nattsericaJ 
weather prediction derived re&vcixvity field, to the 
solution of a system of ordinary differential equations 
sultabte for solving by nrnncrical imttaods. 

The Mcuntcy of tropofiphertc comciioos used by b 
satelUre navjgauon user is 4^}»dt€\i upoA Kbe dcctira^y 
of the numezicn^ weadter prediction (ctmently QA% 
scmdard deviation Car die UKMO mesoscale imdelj a^d 



PO the error intrndueed due to process of correction 
disscmBn&bO^L 

pgr the data set used in this Study, tmn due to tbe 
Oix^^on of a^iTnuthal asytbmctry effects were found to 
be $mall in comparison with micDcrical weather 
prediction errors. Only for usm requiring the highest 
degree of accuracy would aaiiuutlul asymmetry 
xnfbnmtioa be of value. The intact of reispora] 
lesohtioa is significant in relatioa to ibc fiuaeat 
accuracy of the NWp»dcnvcd troposphcric correetioos. 
A ^i^iltcanl additional ciTor is therefore produced rf 
remporal resoludon exceeds 1 hoUTr Spivna] resolutionj^ 
larger ih^n the resolution of the n^csoscalc model, 12 knv 
Are thercfisre found 10 result in signjlicam additional 
ttoposphetic errors. 

For reasons of bandwidth e/^ciency, zenith delay aisd 
n^appine ftuu:tiDn paiaxxAters can he esiprcssed a? 
conections to a pdori &u>del5, 8 Bits coo used to 
dcsctfte irppasphcric zenith d&lay with numnud 
quandsaxioDQt jgioise; a fiirdier 16 B!t9 t»ay be required to 
disseimaaie the elevation mapping ftiaeticza exp^ssed as 
ountinued dOA correction paranteiers. 

The troposphcric correclJoxi maps can be cos^irc^iiid 
with minimal reduction in ixop^ge performance^ but a 
substantial decrease in data size. A coiiabiaacioa of 
efficient quantisation and compression techniques lesidls 
Ul tb^ data-rate for troposphcric coneetion disseDoination 
bvfng suiiafair ibf itaflsmissioa via Global Naviagtion 
SateQiic Systems. 

The zenith delay and mapping fiiftcdon paiamciers can 
be used in asslstuig tbe process satelliiB 
detej ijiuution. 
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P CLAIMS 

L A method of obtaining data &r use in a satclUte positioning system or GN5S 
conqyrisingfhe steps of generating a fest set of less accurate tropospheric delay vahies 
from a first model wiich is known per se, the first model being based on a first set of 
5 parameters, generating a second set of more accurate tropospberic delay values fixnn a 
second model based on meteorological infoimation, and developing a modified said 
set of parameters far us© with said first model so that it then provides a set of delay 
values substantially approximating the second set, and expressing at least one of sadd 
modified set of parameters as a conection paramet^ for communication to a user. 

10 2. A method according to claim 1 wherein said first set comprises tropospheric 
zenith delays. 

3. A method according to claim 2 wherein the first model contama a mapping 
function relating trqpoBpheric delay at a given elevation angle to the zeoifh based Ott 
methodology used by Neil to derive a parametric model of the tropospheric delay to 
15 satellites observed by the user. 

4i A method according to any preceding claim wherein the first modd is based 
on non-meteorological parameters. 

5. A method according to claim 4 wherein said nonrmeteorological parameters 
comprise one^ or any two, or all three of time of year* latitude and altitudo. 

20 6. A method according to any preceding claim wh^ein said second model is 
based on NWP data. 

7, A method according to any preceding claim wherein said second set is derived 
by a ray tracing technique* 

8. A mefliod according to claim 7 wherein parameters generated in the first 
25 model axe generated in the second model from the r^ txacii^ data and the difTermces 

are detemnued Snr a geogi&phloal region. 
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9. A mefljod accoidins to daSm 7 or claiia 8 -wherein a distributiott of 
correctiDii pacametias over at least part or all of the earth»s surfece is derived ft>r 
communicatiaa to at least one geostaticmaiy or nonrgeostattonary satellite. 

10. A method aocotdiag to cMm 9 wherein said distribwtion is subject to data 
redadion to enable a lew c^adty ccmunnnicatiom chamiBl or 

11. A m^od accoiding to claim 10 where said dal^ reduction is performed using 
lossy image compression or sub-sampling techaiques. 

12. A mefiiad accoidiDg to daim 10 or claim 11 and inchiding the step of 
transmitting the reduced data over a ooiiaminioBlioiis channel or date link. 

13. A method according to claim 12 wherean the amomd of data reduction is 
sufficient to pennit said transmissioa within a time substantially lower than (he 
validly time of the said meteorolo^eal inibrmation. 

14. Aineaodaccoidinstoclaiml2whereindieNWPdatodefinesacellsize,and 
if the standaid deviatian for the change in zenith delay for a cell over a period t equals 
the inaccuiBcy in flie said second set, then the data redaction is effective to reduce the 
data to an amount ^sphich csmbe tianamitted over ahandwidfli of 2D0hit per second in 
a time tl substantiany less than £. 

15. A mefhod accoiding to claim 12 wherein V is no more than t/4. 

16. A method of any one of oUdms 12 to 15 wherein said coammmications 
charaiel or data Hnk is to a geostationary or non-geostationary satellite. 

17. A method accoidus to claim 16 vrticidn s£dd gsostafionary or non- 
geostationatysatBlUtorB-traiisniiis at least part of said distribution to a user. 

18. A method accoiding to any preceding claim wherein the or each said 
collection parameter is eixpiessed as the diffetenoe between the parameter of said first 
set and the parameter of said modified set. 
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19. AjjpaiataB of obuuxiiiig data for use in a SEdiettite fiositioning system or C5NSS 
comprising first g^aecatiiig means fisr generating a first of less accurate 
trqpospheiic delay values ftom a &st moSA yrAdck is known per se and which is 
based on non-meteoiological parametew, second genacating means for generating a 
5 second set of more accurate tropoBpherio del^ values flom a second model based on 
meteorological infonnation. and developing means &r develqping a set of modified 
said parameters for use with said first model so that it Ihen provides a set of detoy 
values substantially ^oximating the second set. and said developing means bdng 
Htranged to express the modification to said first parameters as a set of conection 
10 paiametss. 

20. Apparatus according to claim 19 and Biotaddng means fiwconrrrasingsddse 
of correction &ctcas. 

21. Apparatus according to claim 19 or claim 20 and oonqnising transmission 
means for transmitting said set or said compressed set of conectiDn fectots to a user 

15 via a geostetionaxy or non-geostatianary satdlite. 
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